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ABSTRACT 
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A mathematical model of the waterjet propulsion 
plant for surface effect ships consisting of gas turbine, 
reduction gear, axial pump, variable flushed inlet, diffuser 
and fixed nozzle was developed. Control algorithms were 
also developed for the variable geometry flushed inlet and 
for the pump RPM, and were tested by digital simulation of 
the mathematical model. Pump cavitation was included in the 
model and simulation. 


The control scheme for the variable flushed inlet 
involves the maintainance of a diffusion flow pattern about 
the flushed inlet. The control scheme for the pump involves 
the maintainance of pump RPM as dictated by the throttle RPM 
setting. In both control schemes, the rate and hard limits, 
on the hydraulic actuator for the flushed inlet and on the 
fuel flow rate of the gas turbine, were considered. Results 
of real time tests of the developed algorithm are presented. 


Thesis Supervisor: Richard S. Sidell 
Title: Assistant Professor of Mechanical Engineering 


Thesis Advisor: A. Douglas Carmichael 
Title: Professor of Power Engineering 


Thesis Reader: Chryssostomos Chryssostomidis 
Title: Assistant Professor of Naval Architecture 


aa 





TABLE OF CONTENTS 


ABSTRACT 


INTRODUCTION 


I. 


1 ie 


ror. 


WATERJET PROPULSION COMPONENT MATCHING 
I-l General 

I-2 Preliminary Pump Design 
I-3 Reduction Gear Ratio 
MODELING OF S.E.S. AND WATERJET 
II-1l General 

II-2 Fuel 

II-3 Fuel Pump 

II-4 Gas Turbine 

II-5 Reduction Geay 

II-6 Pump 

II-7 Inlets 

II-8 Diffuser 

11-9 Nozzle 


II-10 S.E.S. Mass, Resistance and Waterjet 
Propulsion Thrust 


SIMULATION AND CONTROL OF S.E.S. WATERJET 
PROPULSION PLANT 


III-l General 


a 


14 


15 


21 


Pak 


23 


24 


26 


ao 


41 


49 


54 


56 


Syl 


62 


62 








IITI-2 Dynamic Equations of Waterjet 
Propulsion Plant 
III-3 Pump Cavitation 
III-4 Linearized Dynamic Equation of 
the Waterjet Propulsion 
ZEiI-5 Flushed Inlet Area Control 
ITI-6 Open Loop Simulation to Fuel 
Demand Signal 
III-7 Pump RPM Control 
III-8 Antiwindup Logic for Pump RPM Control 
III-9 Control Algorithm Test Result | 
CONCLUSIONS 
RECOMMENDATIONS 
REFERENCES 
APPENDIX I 


APPENDIX II 


APPENDIX IIT 


64 


72 


74 


78 


82 


89 


2D 


100 


106 





INTRODUCTION 


The inherent geometrical flexibility of Waterjet 
Propulsion system commands a significant interest for its 
use in Surface Effect Ship (S.E.S.). In particular where 
the power level of the propulsion train is high, waterjet 
system, despite its relatively lower propulsive efficiency 
compared to the controllable reversible pitch supercavi- 
tating propellers (Fig. 1) is considered to be a more 
reliable system to fit the sidewall configuration of the 
S.E.S. A great risk is involved when propeller shafts at 
power levels of 20,000 H.P. and above are turned around 
corners to drive propellers. 

Thus the current U.S.A. effort for the projected 2,000 
ton S.E.S. having a cruise speed of approximately 80 knots 
focuses attention to the waterjet system as the best choice 
for propulsion. As a result of the experience of the 100 
ton waterjet driven S.E.S., it was recognized that a vari- 
able inlet geometry was needed to maintain high pressure 
recovery at all speed levels, particularly at hump speed and 
at cruise (Fig. 2). Fig. 3 shows the typical speed resis- 


tance of the S.E.S. One of the design requirements is that 
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at hump speed the waterjet system must be able to deliver 

a maximum thrust of at least 20% greater than hump total 
resistance for an acceptable acceleration performance. 

This implies that the water flow rate for the system at 
cruise speed and hump speed will be approximately the same. 
Since the ratio of cruise speed to hump speed is approxi- 
mately 3, the flush inlet area must have to be adjusted 
accordingly so that a diffusion flow pattern will be main- 
tained always from the region in front of the flusn inlet 
to the pump inlet (Fig. 4). If not, then pressure recovery 
will be seriously degraded and the pump may cavitate due to 
insufficient suction head pressure. 

At present, work is being conducted at M.I.T. to 
optimize the design of components of the waterjet system for 
the S.E.S. so maximum payload can be obtained for a given 
displacement range, cruise speed and hump speed. 

The purpose of this paper then is to model and simu- 
late the dynamic performance of these conponents and to 
identify and bring out problems that may arise in the overall 
control of the system and possible criteria for the selection 
of the components. The case study presented herein is that 
of a 2109 ton S.E.S. whose waterjet system consists of six 


LM2500 gas turbines, reduction gears, axial water pumps, 
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variable inlet geometry flushed inlet and fixed geometry 
nozzle. The simulation program has been set up in such a 
way that thrust reverser characteristics can also be added 


to analyze the backing performance of the craft. 





CHAPTER I 


WATERJET PROPULSION COMPONENT MA'TCHING 
I~l1 GENERAL: 


In order for the waterjet system to function effi- 
ciently, pump design point must match accordingly to the 
thrust requirements both at hump and at cruise speeds. This 
includes the determination of the following: design range 
of variable inlet, area design pump RPM, design pump flow 
rate, design nozzle area. Using the momentum equation the 


thrust delivered by the system is: 


Tg = “g (V, cose - V,) (an) 


where: Ty = thrust develop (lbf) 
P = water density (lbm/f£t?) 
= water flow rate (£t/sec) 
g = gravity (32.2 ft/sec) 
V. = waterjet velocity (ft/sec) 


io= S.E.S. vellocity (f£t/sec) 


Q 
1 


depression angle of the nozzle 
For most design cos@m® is approximately equal to unity 


thus can be neglected in the preliminary calculation. The 





thrust develop then is: 
QV V3 
7. sea & = 1] (222) 
The thrust required at hump speed is: 


Ton — ilo” (Roe + sil (1-3a) 


where: T thrust required at hump 


= speed (lbf) 

R. = bare hull resistance at hump 
Hi speed (lbf) 

Rin = inlet appendage drag at hump 


speed (lbf) 


The factor 1.2 is the required minimum acceleration 
margin. 


fhe thrust required at cruise speed is: 


TR = Rea + Ric (1-3b) 
where: Tao = thrust required at cruise 
speed (lbf) 
Rac = bare hull resistance at cruise 
speed (lbf) 
Rycq = inlet appendage drag at cruise 


speed (lbf) 


The inlet appendage drag for equations (1-3) is: 


R= c, £2 


= <4 
1s! DH 2 Von a 





(1-4b) 


I 


where: C and C inlet drag coefficients 


a ae at hump and cruise 
speeds respectively 
ee and exe = S.E.S. velocity at hump 
and cruise speeds 
respectively 
QO, and Qn = water flow rate at hump 


and cruise speeds 
respectively 


Equating thrust delivered and thrust required and 


substituting equations (1-4) to equations (1-3) results in: 


Oi ; (1-5a) 
H VjH 
g Vou | Vou — (1+0 -6Cpy | 
ony es 
@ thee Vic | 
g Voc | You — (1+0.5C,,) (1-5b) 


The pump head required at hump and cruise speeds from 


Bernoulli equation are respectively: 





2 2 
V V- 2qH 
OH 7H + 4<Q94Iwy, 
H — ae, = 
PH 2g oe) one a a | (1-6a) 
2 2 
H s=} Sener = +—— _ | (1-6b) 
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overall efficiency of inlet 


where: 
"Loa and diffuser 


an 
II 


height of the nozzle above the 
waterline (ft) 


The effective horsepower required (EHP = P HpQ) for 


hump and cruise speeds from equations (1-5) and (1-6) are: 


V5H 2 i 2gHwr 7 m 


ca) = ComSoe V5H ae 
| vee - (1+0.6 ca) | 
Vac\2 2gH 
RaV vs) aig 3 | 
FHP = c70c L(g)" + Roe - oad (1-7b) 


oF Ee | ic ~ (340.5 «| 


Yo 
Note that equations (1-5), (1-6), and (1-7) are all 
expressed as a function of jet velocity ratio (V5/Vo) - AL gh 
the case of a fixed inlet area the jet velocity ratio and 


flow rate at hump and cruise speeds can be related to the 


nozzle jet area (A. ) by: 


A. eoereso 73 
ee atin Vic roe 


Fig. 5 shows a typical propulsive efficiency 
(RV/EHP ) versus jet velocity ratio (V/V) for different 


values of loss coefficient Ky where total head loss 
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A = KV, °/29. Maximum efficiency is observed to occur 
for most practical designs between V5/Vy of i econ? sO . 
Likewise equation (1-7b) can be differentiated with 
respect to (V./Vo) to obtain minimum EHP. However, the 
result may not be compatible to the requirement of equations 
(1-5a), (l1-6a) and (1-7a) because of the fixed nozzle area 
constraint unless the total hump resistance (Roo, + Roa is 
lower than some critical value. For the case where the 
total hump speed resistance is greater than this critical 
value the approach in finding the optimum cross sectional 
area of the nozzle and jet velocity ratios at hump and 
cruise speeds is as follows: 

a. Choose some multiple integral number of existing 
prime movers with a given shaft horsepower 
(SHP) rating. 

b. Multiply the maximum intermittent SHP available 
by the reduction gear efficiency and pump 
efficiency to obtain the maximum intermittent 
total EHP available. 

c. Substitute in equation (1-7a) and solve the 
quadratic equation for (V, vV_.). The two 


OH 


values of (V5 Vow? represent the maximum and 


minimum values of (Vi Vo) for which equation 


(“Jette ibe satistied . 
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ad. Solve for my using equations (l1-5a) and (1-8) 


using values of Vin’ Yo within the range 


H 
obtained in step c. 

e. Solve for Vac/VYoo and check equation (1-7b) 
if the total sustained EHP is within acceptable 
limits. 

£f. Iterate steps a to e until optimum number of 
engines, ee propulsive efficiency is obtained. 

For the case where the total hump speed resistance 


is less than the critical value as dictated by equations 


(1l-7a), (1-7b) and (1-8) the procedure becomes easier since 


Soy d only to be set t ee 
: 2(V3Cc/Voc hee Only O SSS. © Zero O SOLVE Or 


optimum (V5 0/Voc) and the rest follows. 
I-2 PRELIMINARY PUMP DESIGN: 


The procedure stated herein applies primarily to 
determine the rated pump RPM, pump diameter, annulus area 
and other related pump characteristics based on designed flow 
coefficient (%) and head coefficient (), designed flow rate 


(Q) and designed pump head (H,) - By definition: 





= Nain D> (1-9) 








_ SHp _3ep 
Y = U2 ise (1-10) 
60 


where: C.. = water axial velocity (ft/sec) 
U = blade tip velocity (ft/sec) 
Q = water flow rate (ft/sec) 
A, = pump annulus area (£t7) 
N = pump RPM 


D = pump tip diameter (ft) 


The axial pump annulus area (a) is: 


qT 2 
A = = 


2 
= (1 - (4/D)?) (1-11) 





Where: d= pump hub diameter (ft). 


The design flow rate from equations (1-5) and the 
design pump head from equation (1-6) can be substituted 
in equations (1-9), (1-10) and (1-11) to solve for A,, N, D 


for a given value of hub to tip diameter ratio. 


I-3 REDUCTION GEAR RATIO: 


The purpose of the reduction gear is to match the 
Beint of best efficiency of the pump to the point of best 
efficiency of the prime mover. Engine performance maps are 


usually expressed as plots of break horse power versus RPM 
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for various levels of specific fuel consumption (Fig. 6 and 
Fig. 7). Performance maps are also expressed as plots of 
torque versus RPM for various levels of fuel consumption 
(Fig. 8 and Fig. 9). In either case the point of best 
efficiency corresponding to power level required in equations 
(1-7) plus additional power for frictional losses and pump 
efficiency can easily be pinpointed. Thus the reduction 

gear ratio is determined as the ratio of the prime mover 

RPM at this point of best efficiency and the designed pump 


RPM (N). 
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CHAPTER II 


MODELING OF S.E.S. AND WATERJET PROPULSION PLANT 


II-l GENERAL: 


In order to simulate the interaction of the com- 
ponents of the waterjet propulsion plant with the motion 
of the S.E.S., it is necessary to obtain the mathematical 
model of such an installation. The model would then allow 
the development of the appropriate control algorithm with- 
out actually employing the hardware involved. To obtain 
such a model a thermodynamic approach was chosen. Fig. 10 
shows the power flow block diagram of the waterjet pro- 
pulsion train. At this point it is recammended for the 
reader to take a look on ref. (1) and ref. (2) for the 
bond graft technique which will be used frequently in this 


paper. 






FUEL £ _| FUEL GAS CE 
TURBINE | W, 


¢ | FLUSHED INLETI. “pb | wATERJET 
RESISTANCE }| V, AND NOZZLE 0 PUMP 


Figure 10. Block diagram of waterjet propulsion train 





Ze 


where: W._. = fuel flow rate demand (lbs/hr) 
W. = actual fuel ae rate (lbs/hr) 
©. = power turbine torque (ftlb) 

aw 


= power turbine speed (RPM) 


t, = pump torque (ftlbs) 
Ww, = pump speed (RPM) 
a = pump total head (ft) 


t= thrustadis) 


V. = S.E.S. velocity (ft/sec) 


The above diagram is interesting from the point of 
view of thermodynamic accounting of energies involved. The 
light broken arrows represent signal power level and the 
darkened arrows represent high power level transmission. 
The product of the variables affixed to the arrows multiplied 
by some constant represents power. For simplicity, the 
energy storage and energy dissipator elements associated 
with each component and power transmission are not included 
in the detailed modeling of each component which follows. 
Notice that all the energy associated with the fuel will 
in the end be dissipated into the environment primarily due 
to the S.E.S. resistance and secondarily but not avoidably 
due to the "resistances" inherent with the inefficiency of 


each component and inefficiency of power transmission. Thus 





Zo 


from the thermodynamic point of view, for the model to be 
valid in the steady state the power ‘eit for each com- 
ponent must be greater than the power output in either 
direction. Likewise the first law and second law of 
thermodynamics must be satisfied during the transient 


conditions. 
TI-2 FUEL: 


The purpose of fuel is to provide energy primarily 
in the form of heat to the driving engine which is a gas 
turbine in this case study. The energy associated with the 
fuel can be expressed as heating value per unit weight (h). 
Fuel reaction with air mainly is as follows: 


di 


FP 
+ + 
2 CoH + 21 oa 12 CO 18 HO + 6 0, oa 
12 co. + 18 HO + Q) a Q., (2-1) 


where: 1 


FP flush point temperature 


heat evolve in the combustion 
process respectively 


iO 
NO 
II 


The degree on which the above combustion process proceeds 
to the right is dependent on combustion temperature, fuel-air 
ratio and time on which the process is allowed to proceed in 


the combustion chamber of the gas turbine. 





24 


Thus the total heat evolved per unit fuel weight as 
function of time (h(t)) can be assumed in block diagram 


to be: 





= heating values per unit weight 


where: GG = 
of fuel from equation (2-1) 


h(t) = total heating value of fuel per 
weight | 
We = fuel flow rate 
C11 Ts = the exponential time constants of 


the above combustion process. 


cy and T » are very small terms but significant when con- 
vection time from fuel spray area of the combustion chamber 


to the turbine inlet is considered (Fig. 16). 


Ii-3 FUEL PUMP: 


The purpose of the fuel pump is to drive the fuel 
from the fuel tank into the combustion chamber of the gas 
turbine which is at a higher pressure. For the case of a 
positive displacement fuel pump (Fig. 9) represents the 


power flow diagram. 
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Figure 11. Pump power flow diagram 


where: T an = input torgue into the pump 


QO = pump speed 
R,, = pump mechanical frictional resistance 
3 = total lead developed by pump 
We = weight fuel flow rate 
Re = fluid resistance associated with fuel 
value and fuel pipings 
Tw p = pump torque 
I. = fuel inertance associated in the fuel pipings 
m = modulo relating Ge H or modulo relating 
Op 
@W to W 
£ 
TF = ideal transformer of mechanical power to 


fluid power 


1-junction = point of same W or the same We. 
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Since I. and I, are not independent for positive 
displacement pumps, I, can be referred to the fluid side 


(or I, can be referred to the mechanical side). Similarly 


Rn can be referred to the fluid side. For the case where: 
a mW. = eng Val = m 2-2a 
0) rad a a (2~2b) 


p) 
(c) R,= cw. + ——_——"£.__ ) Fes 


£ 2gC. (P Ao) 


then 
92-3) 


2 mikes. 2 1 
(mT, +p) Ge = mig - mo bw, ~ Wel We] (C+ 2§Cz (PA) 2)? 


where: b = mechanical frictional resistance 
coefficient 
7 fluid frictional resistance coefficient 
in pipings 
A, = orifice area of the fuel valve 
Cy = orifice coefficient of the valve 


Ii-4 GAS TURBINE: 


Fig. 12 shows a schematic diagram of a typical air- 
craft gas turbine derivative. It is composed of two major 
units - the gas generator and the power turbine. The gas 


generator is a twin spool engine, consisting of low pressure 
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compressor (L.P.C.), high pressure compressor (H.P.C.), 
combustion section, high pressure menace (L.H.T.) and 
low pressure turbine. The power turbine (P.T.) is adapted 
to the gas generator by an annular diffusing transition 
duct. Thus the operating speed of the power turbine is 
independent of the gas generator speeds. The engine also 
comes with a fuel control package to ensure that starting 
and operating schedules of the gas turbines are met auto- 
matically without the gas generator going into surge and 
without subjecting the turbine inlet to excessive gas 
temperature. 

In order to coordinate the gas turbine package to the 
waterjet propulsion system, a mathematical model is developed 
as a part of this paper. 

For a satisfactory model, it can be assumed that the 
main accumulator of energy occurs in the gas pipe volumes 
and in the rotating parts. Although some heat energy is 
stored in the metal parts, its accumulation and discharge is 
carried out sufficiently slower. Thus during sharp load 
changes the transient is mainly from the dynamic interaction 
of rotors, volume and power turbine speed regulators. 

When a sharp load torque decrease occurs, the power 


turbine speed regulator initial action is to decrease quickly 
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the fuel flow to low value open up discharge value to 
release air intensively from gas en cree to evacuate 
the large quantity of stored energy. After which the 
discharge valve is closed and fuel flow is increased to 
the desired RPM. 

When an increase load is desired on the power turbine, 
it is necessary to increase the gas generator speed and 
therefore the stored energy therein. However in this 
transient process the compressor speed lags behind the 
pressure and temperature increase due to fuel flow increase. 
Thus not only for the protection of turbine against ex- 
cessive temperature but also preventing compressor insurge, 
the rate of fuel flow increase is strictly limited in 
accordance with air flow increase. Surge is also exciuded 
at low speeds by discharging air at the first group of stages 
or from compressor delivery valve. 

In order to see the baSic principle underlying surge, 
consider a multistage axial compressor having the same flow 
coefficient @ = C/U and air axial velocity Cc for all 
stages, at design point. At very low speed Cy. has to increase 
greatly through the compressor since the outlet air density 
1s now approximately equal to inlet air density. This means 


that the first stages are operating at a very high incidence 
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angle or positive stall value while the last stages will 

be at negative incidence angle or Ree aeg en stall value or 
turbining value. Now, with sudden increase of fuel flow 

the compressor outlet pressure will also increase very 
quickly. Thus - also decreases shifting the positive stall 
region toward the last stages. Since the first stages are 
not contributing any pressure increase, back air flow 
rapidly occurs as soon as stall hits the last stage with an 
accompanying explosive sound. 

With the above qualitative background, Fig. 13 is 
shown as the bond graph model of the gas turbine unit. 
Three variables are needed to describe power flow of com- 
pressible gas: mass flow rate (mo), total pressure Coe 


total temperature (T). The subscripts are evident from 


Fig. 12. The rest of the notations are: 


XG, 


torque 

©) = angular speed 

We = fuel weight flow rate 

R = frictional resistance for shafts or orifices 
I = polar moment of inertia 

H = total head of fuel 


V = effective air volume capacitance 
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From the diagram the following are assumed: 


ail) page ey, 7 ont ZOO ~ wee! te) evil oA a st 


i 1 = Pp T =f m.=m 
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"so = Foie weer Ueno Sate Lael 


Above assumption implies that pressure drops along ducts 
and diffuser sections of the gas turbines must be included 
in the overall efficiency of the adjacent compressors and 
turbines. 


From Newton's Law: 


diy, 


Image = (U1 - & LL - by) (2-4a) 


dW) 1 





da@ 


eee Cn = Ca 6G) .) (2-4¢) 


PT PP Pp. P 


Assuming air follows ideal gas law: 


1 RT (2-5) 
then: 
dP 52 . . . 
EE ana — a — 
Ny at (m_. m3) RT. m RT, (2-5a) 
dm, 


a5 = Mo _ m3 (2-5b) 








dP oa ; ; ; ' 
V. es (m, - mix Tey RT + Ley 
where: 
dP 
V5 “der = (hj, - the) AT + mai 
where: 


Mas ~ Ca PA, a\ 05 - Po6)29 


dm3 : 
dt ~ “a5  °a6 
dP 47 s ; ° 
AE en ee aye ag) ee Gee 
dm, = . 
at ~ “a7 ~ “as 
aP 410 : ° 
6 at | (Mag M10) QTD 10N “Se old 
am 
gt aS Calo 


(2-Sc) 


(2-5a) 


(2-Se) 


(2-Sf) 


(2-5g) 


(2—Sh) 


(2-5i) 


eo) 


Equations (2-4) and (2-5) are the equations governing 


the dynamics of the gas turbine in Fig. 12. 


From compressor and turbine maps Poo’? 


Versus 


(mai To1)/(AyPo)) for different values of =— and polytropic 


ol 
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efficiency are plotted (Fig. 14 and Fig. 15) Thus 
ms; can be found as functions Poo: Pol: Je and annulus 
area (A,)- 

From Euler's equation, the torque on the rotor of 


an ideal lossless compressor or turbine is 


7 mr ( &Ce) 
g 
In terms of power: 
— mo rAC, = me Aho 
Therefore: 
v : Moi Ano = 


Equation (2-6) will be true for an adiabatic compressor 
having no storage elements. Thus for a compressor having 
a polytropic efficiency 


( v - 1) A 


@polyc _(n - 1) (n - 1) 
Nn n 
therefore: 
(n - 1) R 
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where 
(n-1) 
Po2 1 T62 
Pol —S eal ai) 


T 
o2 
Ah, = Cp AT, = C61] 3 - 1} - cya |( 222 


7 502) “Prete! ap oid, 
no em mle . 


From equation (2-6) 


v a CT mo? 
Lc 0 "avp'ol Po] 


R 
(C5 ™polyc) | 
) eel (2-8) 


Likewise for adiabatic turbines having no storage 


elements and with a given M polym’ 


(n-1) 
_ —— {jay 
poly: ee as n S 





poly,, 


therefore 


—-Trp = MCT 


a8-p 08 (2-9) 


Mlpoly 
(ey 4 


Equations similar to (2-8) and (2-9) can likewise be 


developed for the H.P.cC., H.P.T. and P.T., thus defining 


the differential equations (2-4). Since from equation (2-7) 


the total temperatures and from performance maps mass flow 
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rates can be found and the rate of temperature changes can 
be reconstructed from the solution of the differential 
equations, equations (2-5) are then defined if the increase 
in total temperature inside the combustion chamber due to 
fuel energy is known. 

Fig. 16 shows the schematic diagram of the combustion 


chamber. 


fp.c. “a5 —- — = =m AG MMP. T. 
a SNS eas YY A hi 
eae 
WwW 
fin 
lof 


Figure 16. Schematic diagram of the combustion chamber 


From the first law of thermodynamics: 


mc T + wfc T + h(t)] = AR 0 
a5 Pq oO5 ‘mas eae ( | ( 
Cc W T + C W - W )T + cC (m_-m _)T 
Dp nout 1OO Peg. fiN fou (oo p, a6 as! 06 

where: C, = specific heat of air 
a 
oe = specific heat of liquid fuel 
it 
C = specific heat of burned fuel 
Peg 
h(t) = instantaneous fuel heating value as 


defined in equation (2-2) 
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Since an and T. in equation (2-2) are very small 
quantities, it can be assumed that h(t) s Whe where hy 
is the heating value of fuel per unit weight, since com- 
bustion chambers are usually designed to carry out complete 


combustion way up before turbine inlet. Thus, equation (2-10) 


mn terms of fuel air ration W/m, reduces to: 


5 
We 
“p_to5 + AS lc, uae + he| 
756 = as. (2-11) 
hee Cc an Be: 
a5 Peg Pa 


TO6 is a significant term since this is the highest temper- 
ature which the turbine inlet is subjected to. For control 
purposes, the measured 756 should account for: 
- Cas | 

1. transport delay , e , where s is the 
Laplacian operator and t 5 is the time delay for the 
burning fuel to reach the turbine inlet, and 

2. time delay due to the inertia of thermocouples. 

It is worthwhile to mention here that among the 
quick acting means of temperature protection, fuel-air 
ratio limiters are widely used which are especially effective 
at starting and sharp load increase. However, fuel-air 
ratio limiters have low accuracy since changes in specific 


heat of fuel or changes in fuel gas ambient temperature and 


pressure do not make setting corrections. Thus there is 
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still a need for very quick acting thermocouples since the 

time constant for heating of the sharp trailing edge of 

turbine blades is in the order of one to three seconds. 
For the case of the LM2500 gas turbine, which has 

a single spool gas generator, equations (2-5a), (2-5b), 

(2-5g) and (2-5h) are simply deleted in the mathematical 


model. 
II-5 REDUCTION GEAR: 


As mentioned in Chapter I, reduction gear is the 
impedance matching component to assure that the power 
turbine rotor of the gas turbine and the pump rotor 
operate with minimum fuel consumption at a steady state 
on power levels most often used. Thus it can be represented 


in bond graph model as in Fig. 17. 


vi 


« 


T- p 
POWER TURBINE pee TR —_——j4 | -————————">—s PUMP 
Wy oe x We 


Figure 17. Reduction gear bond graph model 
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where: TC mp = power turbine torque 
Wp = angular speed of oouee turbine 
Tp = pump torque 
a. = angular speed of pump 


m= modulo of the transformer TF or the 
reduction gear ratio 


I = combined polar moment of inertia of 
pump and power turbine referred to 
the pump shaft 


Re = shaft bearing frictional resistance 
R = reduction gear frictional resistance 
RG ; : 
associated with reduction gear 
efficiency 


Assuming the reduction gear efficiency and shaft 


frictional resistance to be of the form: 


-a) 
fend od 
al 
—— + = —_ 
Were Ones O22 (10 e ) (2-12a) 
—_ “ti = —_ 
Re Co C., W (2-12-b) 


and the constants Ci, C, and C2 are to be determined from 
data, then from Newton's Law: 


I sane = Treen es SRE Te (2-13a) 


p 
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For the case where the power turbine is developing 
negative torque at idle, fuel flow rate during transients 
equation (2-13a) should be: 


m Trp 
I “dt = Wr - Re - C. (2-13b) 





Ii-6 PUMP: 


Pumps are energy converters transforming mechanical 
power to fluid power. In order to understand the basic 
mechanism of pumps, consider the velocity diagram repre- 
sented in Fig. 18 for axial and/or centrifugal pumps 
without guide vanes. The solid lines represent the design 
point velocity diagram while the dotted lines represent the 


off design velocity diagram. 


U = rotor blade velocity 


Wi = relative velocity of entrant fluid with 
respect to rotor 


Cc, = absolute velocity of entrant fluid to rotor 
and also exit fluid from stator 
Cc. = absolute velocity of exit fluid from rotor 


or entrant to stator 


W, = relative velocity of exit fluid from rotor 





C 


rweererwewee weet weet we SF we WS we Se a= = 





STATOR ROTOR 
VANE VANE 


Figure 18. Pump velocity diagram 
For a fixed geometry pump, referring to Fig. 18: 


OC, W, (cos B - sin B tan? ) 


II 
ll 


Q(B ) (2-14a) 


ACY = Wy (cos g'- sin A'tan oo) = 2 (g') (2-14b) 
where: 
c 
W. = v2 +c,7 sin 8 = ———— 
1 1 
1G See 
Ii 
Ww! a= ne + [ey 2 . | C7 
i ea ik sin fo) = 
Ce", 


From the above equation it can be seen that the 
change in the tangential velocity component ( 4C,) of the 
fluid at any operating point is a strong function of 


dynamic similarity "flow coefficient" angle B and a weak 
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function of angle 8 since the deviation angle of the fluid 
at blade trailing edges is approximately constant. Since 
sin B and cosf are periodic functions, then pump and/or 
turbine characteristics must necessarily also be periodic 
functions of dynamic similarity "flow coefficient" angle # . 
From Euler's equation, torque (@ ) and pressure (P) 


ares 


T= poarac, where r is the pump radius 


P= pr Aaciw where W is angular speed 


Substituting equation (2-14) results in: 


I 


v =[p(p)] - 9=-®(B) - 0 (2-15a) 


= | omens) || oe ah ( ) 3 Ce (2-15b) 


Thus, unlike positive displacement pumps and turbines, which 
are transformerly coupled, centrifugal and axial pumps and 
turbines are gyratorly coupled by Rg (6 ire 

If appropriate power dissipation terms, P,, are 


added to equations (2-15) such that: 


= = 0 a 
Py Ty P82 Zz. (2-16) 


where: subscript = actual pump torque and pressure 


Pp 


Pq = power dissipation 
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« 4, 


then for a physical pump we have: 


II 


ja Bg BO Ba By ° Se ey ce (2-17a) 


II 


= " @- (GOP = ‘ _ 
Po ae Ps Rg Ro Q (2-—17b) 


dissipation torque 


II 


where: Va 


= dissipation pressure 


J 
oF 
I 


The bond graph from representation cf equation (17) 


is: 





Ram Lom Ty a 
aC . Pa 
mcommm come oe on 
Ra 
| B1Q 
Ru Re 


Figure 19. Axial or centrifugal pump bond graph model 


polar moment inertia of rotors 


II 


where: Tn 


fluid inertance in pipings 


eH 
Fh 
I 


bearing frictional resistance 


ae 


Re pipe frictional resistance 


II 


Note that the bearing frictional torque a is not the 


Seeeas the pump torsional resistance R,,. By referring to 
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the typical centrifugal pump performance curves of Fig. 20 


E50 


ae” 


g|© 





Figure 20. Typical centrifugal pump characteristics 


and taking ncte of the points (x) where the power flow TW 
or PQ are zero across the gyrator (GY), it is evident 

. 2 
that the dissipation torque Tg = (Geer a) ) and the 
dissipation pressure Pq = P4(Q, W*). It can also be shown 
Enat R,, and a are not independent and can be related by 
some function R,(Q,W). There the bond graph of Fig. 19 


meduces to Fig. 21. 
r | 


Pees 


“| 
| ac 


Figure 21. Axial or centrifugal pump model 
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nae 


i 


where: R R,(Q,0 ) = two port pump resistance 


bearing frictional resistance 


? 


ve) 
Il 


£ frictional resistance on pipings 


From Newton's law and the above diagram: 


In at = sy - Te. - (R,2 “+ R,,@) (2-18) 
qo _ = = - 
Ie ae = (Rw RoQ) P ae (2-19) 
where: % = pump input torque 
€- = shaft bearing frictional torque 
P. = piping resistance friction pressure 
BD = pump pressure load 


ror S.E.S. waterjet propulsion, the axial pump 
characteristic used in the system study (Fig. 22) was taken 
from ref. (6), where the head coefficient YW = glp/U~ is 
expressed as a function of flow coefficient Bp = C/U. Here 


the pump pressure a of equation (2-17b) is: 





2 pe 
bie NY 
H, = U¥ or P = 7" =R_@W - RQ (2-20) 
g Pp g J Q 
where: p = density of fluid (1b/ft>) 


gravity constant (32.2 ey aeeo) 


te} 
II 
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The design @ and Y of the above pump are 0.108 and 
0.51, respectively. To obtain Ty of equation (2-17a), a 


dissipation function coefficient equivalent to equation (2-16) 


is assumed to be of the forn: 


2 
= + ~ a 
min 
where: P, = dimensionless power dissipation 
Po . = minimum Ps (at design point) 
min 


?5 — row cOChotcdoentsat ClSign point 


C, = some constant to be determined 
relating to the efficiency charac- 
teristic of the pump 


Combining equations (2-21), (2-16) and (2-20), then 


Tp and P, of equations (2-17) in terms of pump dimension 


Pp 


reduces to: 


a PAn[ Ng]? 23 
T, = (¥9+ Pg) Geo" D (2-22a) 
No 2 
ev [ov 7 >| 
ee ae (2-22b) 
P g 
where: oe = pump torque (ft-lb) 
Y = dimensionless head coefficient 
Pp = dimensionless flow coefficient 
Ps = dimensionless power dissipation of 


equation (2-21) 


P = water density (64 lb/ft?) 





49 


a 4, 


A, = pump annulus area (£t*) 
N = pump RPM 
D = pump blade tip to tip diameter (ft) 


II-7 INLETS: 


Flushed and semiflushed inlets shown in Fig. 23 
and Fig. 24, respectively are considered to be more 
practicable to use for S.E.S. than ram inlets. 

The geometric parameters affecting the design 
pressure recovery and drag characteristics of flushed and 
semiflushed inlets are: entrance width to depth ratio, 
ramp wall divergence and ramp angle. For fixed geometry 
inlets, performance at various operating points is best 
described by efficiency of pressure recovery Ms. asa 
function of velocity ratio Vs/V, where ve is the velocity 
of the craft and V; is the average axial velocity of the 
fluid across the inlet throat area A;. Some typical pressure 
recovery performances of flushed and semiflushed inlets 
(ref.(32)) are shown in Fig. 25 and Fig. 26. 


Similarly, the drag coefficient 
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FIGURE 23 FLUSHED INLET 






STATION 2 
STATION 1 


FIGURE 24 SEMIFLUSHED INLET 
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FIGURE 24 FLUSHED INLET DRAG 
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can also be expressed as a function of the ratio Vie? Vn 
to describe the drag (D) ee eet eee the inlet 
at various operating conditions. Some typical drag 
characteristics of flushed and semiflushed inlets (ref. (32)) 
are shown in Fig. 27 and Fig. 28. 

From the above pressure recovery characteristics 


the bond graph representation of inlet is: 


1, 


wale 


2 
Pe een ey eed ee 
at 24 q's) @ 


| 


Re 


Figure 29. Flushed or semiflushed inlet bond graph model 


whe re: I. = fluid inertance in the inlet section 
Pol = total pressure of the fluid way 
ahead of the craft 
Re = equivalent resistance of the inlet 
plus the boundary layer losses at 
the sidewall of the S.E.S. in front 
of the inlet 
P 59 = total pressure at the exit of the 


inlet or at the entrance of the 
diffuser section so that at the 


tead tate: P = .P 
Bee So OZ aL oul 
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II-8 DIFFUSER: 


The purpose of the diffuser is to increase the 
static pressure at pump inlet by converting the kinetic 
energy to static head. For a given expansion angle and 
fixed inlet condition, the total pressure loss of a 
diffuser is proportional to the theoretical value of the 
total pressure loss for a sudden expansion of the same 


: ; : Aout 
exit area to inlet area ratio (AL, = 77) 
R Ain 


This pro- 
portionality constant is called diffuser loss factor (Kg) 


and is defined by: 


f Ue sE)) 
AR 
where: Ee = total pressure at entrance to diffuser 
chi) 
P = total pressure at exit to diffuser 
Pout 
a= dynamic pressure at entrance to 
diffuser 
a 
_ (285s 
G1 2g 


Ar = area ratio 


Figure 30 shows the diffuser loss coefficient K as 
a function of area ratio (Ap) and diffusion angle (290) for 


various geometrical shapes of diffusers. The bond graph 


representation, then, of the diffuser is: 
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Figure 31. Diffuser bond graph model 


so that at the steady state: 


II-9 NOZZLE: 


(2-24) 


The purpose of the nozzle is to convert the excess 


static head of the total head by the pump into kinetic 


energy head so that thrust can be developed. 


nozzle, unlike diffusers, has negligible loss so that 


efficiency of 96-99% can be easily attained. 


graph representation of nozzle is: 


A well-designed 


The bond 





oe 
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ev 
Poin — Pa 2g 
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Figure 32. Nozzle bond graph model 


so that at the steady state: 


ev," 
Riga 
2g = %,'Po. - P,) (2-25) 


where P.. total pressure at the pump exit 


in 


rE 
a 


Un 


I 


atmospheric pressure 


il 


nozzle efficiency 


II-10 S.E.S., MASS, RESISTANCE AND WATERJET PROPULSION 
TRV + 
The S.E.S. bare hull resistance data used in this 
study are tabulated in Appendix I. Note that if the 
boundary layer losses due to S.E.S. sidewall frictional 
resistance are included in the conputation for the total 
pressure recovered at the exit of the flushed or semiflushed 


inlet, then by energy conservation, an equivalent resistance 
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must be subtracted from the bare hull resistance data above. 
Otherwise, the propulsion system will een eeeG twice. 
The assumed added mass for the forward motion of the S.E.S. 
was 10% of the total displacement. 

From the conservation of momentum, the instantaneous 


thrust developed for a fixed reference frame is: 


T = =a if [4 vav |+ J [ eae as) | (2-26) 


or in terms of component thrust in x and y direction: 


pl dQ Pe 
joe = +— (V. cosvc - V ) (2-29a) 
I O 


* 





x g dt g 


pt dQ P2 
ju at + cae (Va sin & ) (2-29b) 


4 


| 


where: Li, and L,, are the distances between inlet 


entrance and nozzle exit in x and y directions, respectively. 
At the steady state, the “power input" to develop 


thrust is: 





y) 2 
rF E 
AP, Q= - 1 Q (2-28a) 
in 2g 52 


and power output is: 





i Vo (2-28b) 
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where: 


APQ - TaV. > 0 


d 


Since 


then it is evident that the power conversion of eguations 


eo) 
(2-8) is gyratorly coupled with modulo m, = “~Q7 v, 7 il 


and the bond graph representation is: 


I §.E.9. 


APin 4 LD Fese GY Tg 1 
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0%% 
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Figure 33. Bond graph model of thrust develop in S.E.S. 
waterjet system 


The equivalent dissipator resistance 


AP AP. - AP 
Peeetoss: _ ___in 22) eff 
qd” Q 7 Q 
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where 
ae el . 
AP OFFt g Vee V +m, A 
Pp 2 v.7 
Pp, =+tfy, [te - 2 
aH gl UE: 
© 
is: 
pP 1 W5l V3 ve 
—— as poe ak Ay) ls = 
Ry = "gG |VYol Vo 3 je eo, al7 © 2a) 


The overall waterjet system bond graph model which 


combines all the different components is shown in Fig. 34. 
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CHAPTER L121 
SIMULATION AND CONTROL OF S.E.S. WATERJET 


PROPULSION PLANT 
III-l1 GENERAL: 


For the purpose of controlling the power turbine 
speed, the implementation of the detailed model of the gas 
turbine was found to be unnecessary. This.was primarily 
due to the higher frequency response of the combustion 
process and pressures in the storage volume of the gas 
turbine in equations (2-5) compared to the frequency response 
of power turbine speed. However, for a very low value of 
polar moment of inertia of the power turbine drive train as 
referred to the pump shaft, the poles and zeroes of 
equations (2-4) and (2-5) will be significant in the design 
of a stable high frequency controller. Likewise, if the 
primary concern of designing a controller is the avoidance 
of surge in the compressor and of overtemperature in the 
turbine inlet of the gas generator then the detailed model 
will be necessary. 

In the case study it was assumed that the LM2500 


gas turbine used in the simulation was adequately protected 
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from overtemperature and surge by its own control package. 
This gas generator protective control setup was simulated 
by putting a scheduled limiter on the rate on which fuel 
flow demand signal can be increased. This rate can be set 
to have lower value than the built in limiter in LM2500. 
The purpose of this limiter is to allow time for the gas 
generator rotor to build up speed and in turn air mass 
flow rate in order to neutralize the increase in pressure 
for surge avoidance and the increase in temperature due 

to sudden fuel flow rate increase. 

As for the rest of the gas turbine dynamics, it was 
noted that the instantaneous torque developed at the power 
turbine is a function of the instantaneous air mass flow 
rate and power turbine speed. Air mass flow rate in turn 
is dependent on the instantaneous temperature and pressure 
at the storage volume of the duct between power turbine and 
the last stage of the gas generator. Thus the transfer 
function of the power turbine torque from equations (2-4), 


(2-5), (2-7) and (2-9) can be approximated by: 





Cr 1 i “|= it — = 
Tss | Teagstl)| | Tst1 (Ge Ce + (Cagtts) sti (le 


NX 
tl 


where: instantaneous power turbine torque 


I 


Cog 


steady state power turbine torque which 
is a function of fuel flow rate (We) 
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and power turbine speed (N,,) 
Cac = time constant for the gas generator 


(a time constant for the volume between 
power turbine and gas generator 


s = Laplacian operator 


Since Tae and Ty are in the order of 107+ second, 
which is very much lower than the time constant of the 
power turbine drive train speed, the first term of equation 
(3-1) can be neglected so that: 
(s: 1 
Tss = (Tog + Ty) stl (3-2) 
The plot of age as function of fuel flow rate Oe 


and the power turbine RPM is shown in Fig. 8. 
III-2 DYNAMIC EQUATIONS OF THE WATERJET PROPULSION PLANT: 


With the above approximation, the water propulsion 
system was reduced to the following non-linear ordinary 


differential equations: 


a. The S.E.S. velocity Equation 


Se Cae [,, (Vo) OY) ie Ry (V5/ 29g) | (3-3a) 
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where: V. = S.E.S. velocity (ft/sec) 
W = S.E.S. mass plus 10% added mass (lbm) 
GM = Gravity (32.2 ft/sec’) 


T = waterjet propulsion thrust (lbf) 


x 
R. = S.E.S. bare hull resistance (l1bf) 
= flushed inlet appendage drag (lbf) 
2 
a flushed inlet throat area (ft ) 


The instantaneous thrust (TL) in the longitudinal 


direction of the craft from equation (2-27a) is: 


pl, dQ dv P2 
TT = — Cram - ——— Vy, cosio = V } (3-3a-1) 
x g dt dt g O 


where: P = water density (lbm/f£t?) 


L. = longitudinal distance between the flushed 
inlet entrance and nozzle area (ft) 


A, = effective cross sectional area in 
longitudinal direction (ft*) 


V. = waterjet velocity at the nozzle exit 
(ft/sec) 


O& = depression angle of the nozzle 


The S.E.S. bare hull resistance (R,) should be a 
function of S.E.S. velocity (V, and S.E.S. operating draft 
Or air cushion pressure. In the simulation, however, it 
was assumed that the S.E.S. was at designed cushion pressure 


above hump speed. Below hump speed a resistance (Ro) 
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2 
proportional to uC was assumed, so that: 


Fe= 167.72 V,* 0<V <50.7 (3-3a-2) 


2 
S O O Or 


The flushed inlet drag R from equation (2-23) is: 


Ve 


vo 


BY dt 
ap V 


(3-3a-3) 
O 

For the purpose of controller design the flush inlet 
drag coefficient Cg was assumed to have a sharper point of 
minimum value at V;/Vo=0.-8. The minimum value of Cg was 


also assumed to be equal to 0.15 consistent with the published 


experimental data, so that: 


Cq = 1.0 - (0.8/0.85) (Vj/V,) 1.0<V;/V£0.8 
a > OFZ Vi/Vo - 0.01 V;/V_7 0-8 


b. The Pump Speed Equation 


GN _ 60 v ; ; ; 
dt ~ Tz | te n RGR) TN) T,,(,0)] (3-3b) 
where: N= pump RPM 
I, = polar moment of inertia of power 


turbine, reduction gear and pump 
rotor referred to pump shaft (ft- jase" ) 
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Ce = effective engine torque at pump shaft 
(ft-lbf) 

Tp = shaft bearing friction torgue (ft-lbf) 

Tr = pump torque (ft-lbf) 


The effective engine torque at pump shaft from 


equation (2-13) is: 





= >O = = 
(= RGR) oe Pc (os (3-3b-1) 
(4 i: v 
= (RGR) <0 
aia " RG E 
where: RGR = reduction gear ratio = 5.55 


MRG = reduction gear efficiency 


7 -N, /895. 
‘MrRG 0.75 + 0.22 (1 -e "5 ), where 


Np = power turbine speed 


The shaft bearing torque (7) is assumed to be: 


Ge = 1000 + 15N (ft-1bf) (3-3b-2) 


The pump torque (Vp) from equation (2-22a) is: 





il 


2 
A N 
(YP@+ PR.) : ae n| D (3-3b-3) 


T. 2g 


where 


= 
n 
o = ND 





68 


A curve is fitted to the pump characteristics of 
Fig. 22 for the head coefficient (~) as a function of 


flow coefficient (Pp) so that: 


Y = -83.3° - 1.9169 + 0.85 $0.04 
= 220° Se2cecget 1.308 0.0 <f<0.065 
Y= 68.217 g* + 10.173 + 0.2069 0 .065< $<0.130 


Y= -318.45f° + 72.7136 - 3.6949 0.130 $£0.152 


The dimensionless power dissipation IP, of equation 
(3-3b-3) was normalized to a pump having an efficiency of 


90% at design point so that equation (2-21) yields: 


P, = 0.00623 + 12.4 (¢ - 0.108)? 


c. Waterjet Flow Rate Equation per Pump 


fe adv ) 
= —_ = 


where: Q = water flow rate fee eee) 
3 
P= water density (lbm/ft ) 
; : : ; 2 5 
I. = fluid inertance in pipe (lb-sec /ft ) 
Hoy = pump net positive suction head (ft) 
H, = waterjet kinetic head at the nozzle (ft) 


He, = head generated due to S.E.S. 
acceleration (ft) 
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« A 


The pump head ee from equation (2-22b) is: 


H, = : lor p| : (3-3c-1) 


The net positive suction head Haw from equation (2-25) 


and from flushed inlet efficiency characteristics is: 





2 
as Es “9 3%p - Hy, (S=se=2) 
where: M3 = flushed inlet efficiency 
Up = diffuser efficiency 
Hor = pump height above waterline 
Mm; = 3 CE OAS 
Mp = 1-Ks OR)? 


The waterjet kinetic head at the nozzle (H) is: 





2 
Ma 
where: M x = nozzle efficiency w 1.0 
V5 = O/C _ Anoz = waterjet velocity (ft/sec) 
C. = nozzle coefficient of area con- 


traction ~ 1.0 


Axoz = cross sectignal area of nozzle 
Ecoat. (pce) 
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The head generated through S.E.S. acceleration (H,) 


is: 


L,{dV, 1 dQ 
Hy = “gl|dt - AL dt (3-3c-4) 
where: L = distance between entrance to flushed 


meet and nozzle Exit in the ve 
Girection (ft) 


d. Power Turbine Torgue Equation 
ae 


dt 





= 0.5 | oe Me Me = | (3-3d) 


0) 
KK 
¢)) 
R 
ll 


instantaneous power turbine torque 
(ft-lbf) 


se 
{ 


= steady state power turbine torque 
as found in the engine performance 
map 


a 
il 


power turbine RPM 


The factor 0.5 is the \ Cae ee) term in equation (3-2). 


e. Fuel Flow Rate Equation 


aw 


“qr = 2.0 (W Ww.) (3-3e) 


fd é£ 


instantaneous fuel flow rate (lbm/hr) 


= 
@ 
Ky 
* © 
=, 
Il 


= fuel flow rate demand signal (lbm/hr) 


=, 
Fh 
O 

| 
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Equation (3-3e) assumes that LM2500 control package 
has a built-in servomechanism for the fuel pump with a 
time constant of 0.5 second and a constant gain of 1.0. The 
rate limit and hard limit were incorporated in the fuel flow 
rate demand signal so that 
Ww. <Ww 


fd ~ Se 


<a WwW 


W <wWw 
min fd— fdnax 


fd 


The minimum fuel flow rate demand signal W was 
£a 
min 
set to idling fuel flow rate of 800 lbm/hr and the maximum 
fuel flow rate demand signal Weg was set to 10,000 lbm/hr. 
max 
The maximum rate of change of the fuel demand signal Weg 


Max 


was set to six seconds from Weg . to Wey : 
min max 


f. The Flushed Inlet Area Servomechanism Equation 


dA. 
—= = 0.5 (A:q - A.) (3-3) 
dt ~° 1d i 
where: A, = instantaneous flushed inlet throat 
area (ft~) 
BL = flushed inlet throat area demand CHESS 
a 


Equation (3-3f) assumes a hydraulic servoactuator for 
the flushed inlet throat area with a time constant of two 


seconds and a constant gain of 1.0. In the simulation, the 
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hydraulic actuator was set to have a hard limit of 7.0 £t* 


2 
maximum, 1.5 £t minimum, and a rate limit of 0.5 eae 


so that 


15 < Te < TO 
IiI-3 PUMP CAVITATION: 


To simulate pump cavitation during transient 
condition at various pump operating points, a dimensionless 
suction head criteron was assumed to be: 


2 


Y _ “gue _ a7 


SVC 


() 


SVC 


where: Hove critical net positive pump suction 
head (ft) 


I 


pump blade tip velocity (ft/sec) 


G 
iI 


A curve is fitted into the pump cavitation data 


(Fig. 22) so that ee is defined as a function ) . The 


pump is assumed to cavitate if the actual dimensionless 


suction head Pa Bes fell below YW... The degree 
U 2 


of degradation in pump performance was assumed to be of 


the form: 


Y . =V11 - ¢, (Ysyer¥sua) Y sva < are (3-5) 


SVC 
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where: Yo = dimensionless head with cavitation 


Y = dimensionless head without 
cavitation 


C, = some constant to fit the pump 
cavitation data 


Substituting equation (3-5) to equation (3-3c-1), 


the pump head generated with cavitation is: 


Y baee te 2 
H -t) - — s ai p| Ysva <¥sve 


Pp (e=6) 
The dimensionless power dissipated when there is 
cavitation P, is also assumed to be of the form: 
VYsve 7 Ysva 
PR, =P_{1.0+¢ [ Tevet sva'} Vsun Lf Gop 
Cc qa 2 ae SVA <TSVC 
where: Ps = dimensionless power dissipated without 
cavitation as defined in equation (2-21) 
Co = some constant to be determined from 


pump data 


Substituting > and P in equation (3-3b-3), pump 


dc 


torgue with cavitation is: 





t, 7 Ib, ° Pad = E D° ¥ sva< Y svc 


S| 
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IIi-4 LINEARIZED DYNAMIC EQUATION OF THE WATERJET PROPULSION: 


The linear expansion of equation (3-3a) about a 


steady state operating point yields: 














; P g | 7, a T,, OR. 
‘- _— MSG + W “OV Ee) a0 (o-OR) = OV (VV, 4) 
OR; OR: OR 
- “Av “EVo3 - “FQ (0-O,) - Ha, Ai-Ajd! 


The subscript o represents the steady state operating 
point values of the state variables and the partial deri- 
vatives are evaluated at the corresponding steady state 
operating point values. Similar expressions can be derived 


for the rest of the equations of (3-3). At steady state: 


V ae ro) = We = W aac — 
ote) e) oO Eo ES lo 
Defining: 
8) 6 e Oo 
eee = Ve. and Vem One Ora oe 
OO. els and OC BS Q. 


etc. yields the linear form of equations (3-3). 





where: 


11 
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The block diagram 
solid lines in Fig. 35. 
The cuts in 


controller. 


and 2 will be discussed 























[2.24] 
OQ 0Q 
[2] 0 ‘ 
OA; 
|? ss a “E| = 0.5 ror| 2"ss| 
0 Np Oo N Oo Np 
oe | 
| - Op |= -G.5 
27s | 
0 We 
W 
|- <2 =| 20.50 
QW, 
a i] 
= = -0.5 
[oa 
|e | - a 
OWeg 
Ee 
=O. 5 
dAig 


of equation (3-8) is shown as 
The broken lines represent the 


x 


f 


the diagram at points y 


f 


in the succeeding sections. 
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t 
eo a, 


TIIT-5 FLUSHED INLET THROAT AREA CONTROL: 


For the performance of semiflushed or flushed inlet 


to be optimum 1) pressure recovery must be maximum and 


2) drag must be minimum. Fortunately, these two 


criteria 


were not contradictory, as can be seen in the performance 


curves of Fig. 25, 26, 27, 28 of Chapter II. At 
to ship velocity ratio of V/V, = 0-8, the above 
ditions are satisfied. A control algorithm then 
set up in such a way that V./V, (or equivalently 


O/AV. = 0.8) at steady state and at least close 


the inlet 


two con- 


can be 


to 0.8 


during transients. A preliminary proportional and integral 


control for the flushed inlet actuator was set up in the 





form of: 

As = K, — - 0.8) 

1dp FN AiVo 

a = O : ; 

ara at ‘AL, Yate) 
so that: 

= ‘ + 
Aa Aidp AS 


Substituting (3-9) to (3-3f) yields: 


Q 





(3-9a) 


(3-9b) 


(3-9c) 


® 
A, = -0.5Aj + 0.5Aiq, + EoD (AAV, - 0.8) = (3-10) 





ih) 


Taking the linear form of equations (3-9b) and 


(3-10) and including them in equations (3-8) yields: 











Vv, Ai, 9 AYzZ 0 O AY, O Ve O (saab) 
Ny D Fe A, Ao, 9 O 0 N, O 
Q, Az, Ag A353 Oe On 0 ee Q, O 
can Oe aoe Oem Tn, | +] 0 Mea | 
Wee JO © Oy Fe ee Wey Bay 
| A A AA 
Se pom ace SC GG ie 
idy, 471 9 Azz 9 0 Aye 0 idy, O 
where: 
Q 
= a alive [ QO ] 
Aci = 0.5 Koa 2 Vo = -0.5 Koa AVV2 
Ena 
a ae Ciao 2a [ i 
Ga. ae PA Oel = we Panes 
Q 
ae a[x,vo] QO 
A. = Aggt0-5Kpa—54 = -0.5 Kpal 3 ave 
1 1 Oo 
A =0. 
ay 0.5 





> 
II 


Q 
zal ia | ave] 
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1 
an ea ] 
io 


= Q 
A = -K | 
A 2V 
76 I A;“V. 


This PI control is shown in Fig. 35, in which: 





= 2 
By = AiVo 

1 
ae io 


To obtain the steady state values of V ,N,Q, 
elo O O 
CEO We, and ie an open loop response to different 
levels of Weg using conservative values of K and K 
were obtained. Some typical responses are shown in Fig. 36. 
To obtain an optimum Kian and a an open loop 


transfer function ay by cutting the loop at point x 


* 


in Fig. 35 was obtained using MBODE program. Such open 


loop transfer functions can be represented in matrix form: 
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v, Ay, 9 AY, 0 0 0 0 Vv rT 

NL DP ye CO aaa NU 0 

on A31 As, A330 0 oO O Q. 0 

v e.| =| 9 Aga 0 Aggy Ags 09 Of} | TE | +] 0 Y. | 
We 0 0 0 O A. O Of | We, O 
A Ag, 9 Agz 0 0 AgeA. As 0 
Alas. 71 9 A 0 9 Be Of | arg 0 


< 
= 
Tl 
Tah 
° 
© 
° 
© 
Ie 
Loy 
i 


where A,,¢ of equation (3-11) is set to zero. 

The zero location can be adjusted by the ratio of 
Kra/Kpa depending on the location of poles and zeroes of 
the hydraulic servoactuator for the flushed inlet at various 
operating points. The hard limit and rate limit of the 
actuator on the other hand, sets a constraint on the 


maximum value of Kya that can be used, otherwise excessive 
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overshoot occurs. For the case study presented, the hard 
limit and rate limit on the actuator PR cone by 
setting Kia = 0.5 and K, = 5.0 and putting a hard limit 
on the integrator portion of the controller equal to the 
hard limit of the actuator, so that the maximum and 
minimum values of Aid, were 7.0 and 1.5, respectively. 
Thus for the control of the flushed inlet geometry 
some means to measure Q, V, A; must be provided in order 
to keep track of V3/V.- The time delays associated with 
these measurements must also be included for the proper 


location of the zero of the PI controller. 
TII-6 OPEN LOOP SIMULATION TO FUEL DEMAND SIGNAL: 


The fuel demand scheduled limiter used in open loop 
simulation was set to have a maximum rate of change from 


idling Weg to maximum W in six seconds and a hard limit 


fC) 
of idling Wrg = 800 lbm/hr and maximum Wea = 10 ,OCOmLom/hr. 
' | | PPxfdq 2% 
The Hy term in the flow equation (3-3c) and the g | ata dt 


term in the thrust equation (3-3a-1l1) were neglected since 
they were negligibly small compared to the other terms. 
Since S.E.S. under study uses six engines, the thrust 
equation (3-3a-1) and the inlet appendage drag equation 


(3-3a-3) were multiplied by six for the simulation of the 
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straight line motion of the craft. Computer runs were 
made uSing a stepped fuel demand signal (Weg) from each 
engine idling fuel flow rate of 800 lbm/hr to various 
levels of W,, up to maximum Wea of 10,000 lbm/hr per 
engine. The stepped Nee Signal was filtered by the maxi- 
mum rate change of the Scheduled Me limiter with the 
setting discussed above. For a realistic open loop simu- 
lation, the flushed inlet area control loop discussed in 
the preceeding section was closed. Pump cavitation was 
observed at the stepped Wes of 5,000 lbm/hr/engine and 
above. Some responses of pump RPM, pump flow rate (ft?/sec), 
eep.S. velocity (ft/sec) and flushed inlet area UU) eS: 
the stepped W-, made are plotted Pe 1g. Sr 

Fig. 36 responses, however, were obtained using a 
higher value of fluid inertance (I-) in the pipings than 
the computed value for a faster YSYS solution (1i.e., using 
higher value of time step). Comparison of the response with 
the actual value of I (1i.e., using very small value of 
time step) yields a similar result. The higher value of 
I, used served as a filter to the oscillatory solution of 
YSYS to pump flow rate during transient and especially in 
the cavitation regime when the time step is greater than 


eel second. 
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III-7 PUMP RPM CONTROL: 


Since pump shaft RPM appears to be the easiest to 
measure accurately and appears to be the most effective 
indicator of power turbine stress overload, it was decided 
to base the fuel demand signal on the error (E.) between 
the throttle RPM setting (N,) and the actual pump shaft 
RPM (N). 

By inspection of Fig. 35, the open loop transfer 
function N,/Weg, yields a type O system. Such open loop 
transfer function, represented in matrix form of equation 
(3-11) with the output matrix [c] = [0100000] , was 
verified using MBODE program and the result yields no pole 
at the origin in the s Pai. Thus, to reduce the error 
(Ey, = Np, 7 N,) to zero at steady state, an integral 
control is needed. A proportional control is also needed 
to bring N, close to Np, faster for large error (Ey) - In 
order to obtain the gains Koy Kon for above PI control, 
MBODE program was used to locate the open loop poles and 
for various levels of W i es = 


eae FQ’ 
800 lbm/hr to 10,000 lbm/hr. 


zeroes of N,/ 


In the MBODE solution, however, of equation (3-11) 


the A332 term was so big compared to the Aly Ayo: Anat bee 


and Age terms that a big numerical error resulted. It 
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was therefore necessary to reduce the order of the system 


to six by setting: 


This is equivalent to setting the fluid inertance to zero 


or shifting the pole due A to -ee in the s-plane. Such 


33 
a procedure will yield a solution very close to the true 
solution since the pole due to Bees being so far away to 


the left of the other system poles and zeroes, is effect- 


ively at -c©®. Solving for Q yields: 
* - 


A 
go, = stv, - 2 N , (3-13) 
A A oe 
Bo 33 
Substituting equation (3-13) into V ; ; A. and A. 
- Ey ay ids 
of equation (3-11) results in: 
Vy ayy 255 0 0 te 0 Vy 0 
N. ayy ao5 ae O O O Ny O 
UT 0 a A332 a o) 0 ao 0 
ee! 32 34 E,| Yea 
on O O O ang O O We by 
i. si, 5 (UR Se cy | | car : 
id, ae, 362 0 O age a *idy ail 





where: 


475 


2211 


a3 2 


I 


91 


opt 


yy 


22 A33 


24 


42 


44 


45 


Ag3 431 
61 aie 
R63 A32 
a3 





22 


5 
e @ 


With the above PI control, since W = 2.0(Weg-W,) 
he 


and Wg = (K, 5 7 


transfer function N /N 
* Re 


Vv, ana 
N, sal 
Te 0 
ie. AA hae 
A, 451 
Asay “61 
Mea. 0 


K 





Ss 


N 
) (NL-N), then the closed loop 


is represented by: 


aio O ajg O 0 V, 0 (3-15) 
Qn5 423 O O O O N, O 

Az 833 834 0 O O ee = 

ayo 2a O O ang We + D, *s] 
aso a55 856 0 A O 

Aes 0 Aes 0 0 a 0 

as 5 0 O O O Weg bo 


with the output matrix [c] =[o 10000 0] 


where: 


-2.0 Kon 


= +2.0 


IN 


4+#2.0 K 
PN 


+Kin 








“Ge = Byes 
A 

| “13 31 

61 71 ie 
- Aj3 A30 r 

62 2G * 
“65 ~ "76 
Be 

4 ~ F- 


An MBODE result for the value of matrix A of 


equation (3-14) at Weg 


[c] = [0 1000 0] indicates that three of the zeroes 


= 5,000 lbm/hr with output matrix 


cancel out three of the poles of the open loop transfer 


N, (s)/W (s) and is plotted in Fig. 37. 
k 


qd 


jo 





Figure 37. Open loop poles and zeroes of N,/ Weg 
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Thus in the PI controller of the form: 


: Kin Kpy | San 
ea = + 


a ss Kon 


Kry/K,y can be adjusted so that the zero of the controller 


is in the vicinity of -0.234 and the gain oe also adjusted 


so that closed loop pole is within the 45° line indicated 


in Fig. 38 for a fast response with less than 5% overshoot. 


‘N 
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Figure 38. 


fd 


Root locus plot for K Ko aoe 24 
P In’ PN 
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III-8 ANTIWINDUP LOGIC FOR PUMP RPM CONTROL: 


In order to compensate for 1) the excessive over- 
shoot and oscillation during transient caused by the hard 
limit and rate limit on Weg Signal and 2) the shifting of 
the two rightmost uncancelled poles of equation (3-14) 
towards the origin at low power level, the K and Koy 
were adjusted accordingly in the form as shown in Fig. 39. 

The idea behind such a scheduled Ky and Kx 1s to 
set the gains of the PI control for small error signals to 


be equal to the optimum values as found in the linear 


analysis of section III-7, so that: 


f_ (N_) 
*Pmax 1 R 


K 


Inax 


eo 


The £, (N,) and £. (Np) are the plots of the optimum ee 

and KON as function of N. which will yield a fast stable 
response with less than 5% overshoot about an operating 
point a8 for small error disturbance. For large error 
disturbance, however, the value of ae is too large to 
be compatible with the rate and hard limit of the scheduled 


fuel demand limiter. Large overshoot occurs because by the 


time the integral part of the fuel demand signal (Weg ) has 
I 
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Koy = KomMAX (1-0 - ABS (Ne-N)/100-0) 


Ke min 









ABS ‘4 Ne-N) 


Kin 


Kr MAX 


Krn = Kiwax (1.0 - 489 CNa-N) /100.0) 





ABS (Nea-N) 


FIGURE 39 PUMP RPM PROPORTIONAL AND INTEGRAL 


SCHEDULED GAINS 





oH 


reached the correct value of W for a steady state Nps 


the error signal is still lagging too far behind. This 


is due to the fact that the correcting signal W 


fd acting 


on the fuel pump servo is considerably less than the 
combined integral and proportional fuel demand signal 


+ Weg ) due to the filtering action of the scheduled 


Pp 
Weg limiter. Since for large error the proportional part 


(Ww 
Gl 


of the fuel demand signal (W,, ) by itself is large enough 
P 


to be filtered by the W 


fd limiter, K therefore can be 


IN 
set at a lower value and still get the same speed of 
response in pump RPM (N). Thus oe can be set so that 
by the time the error signal goes to zero, the integral 
part of fuel demand signal (Weg) is equal to the actual 
fuel flow rate (We) for reference pump RPM (N.)- Over- 
shoot through this setup is eliminated. 

However, for K, = constant = K , it was observed 

max 

that some secondary oscillation occurs during the accelera- 
tion transient (Fig. 33b). This is due to the hard limit 
on the Weg signal for large error. When the error signal 
was decreased to a value such that Weg was no longer hard 


limited, the Weg Signal decreased rapidly, which caused 


the actual fuel flow rate We to decrease slightly. To 
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overcome this effect, a negative Sie iva ie Ciaensste)e of Kan 
with respect to absolute value of error was set so as. to 
maintain always an increasing fuel flow rate (W,) during 
acceleration transient. K was also set so that for 


Pmin 
large errors Weg will still be hard limited. 


r 


III-9 CONTROL ALGORITHMS TEST RESULT: 


The following hard maneuvers were simulated to test 
the control algorithms for pump RPM and for the flushed 


inlet throat area: 


(a) Acceleration from 109 to 635 RPM Throttle Setting 


This simulation was conducted by putting a step change 
in RPM throttle setting. The 109 RPM corresponds to idling 
power setting and the 635 RPM corresponds to 95% of the 
maximum intermittent power rating of the engine. The results 
are plotted in real time (seconds) in Fig. 40. The labels 
and scales are: 

V = S.E.S. velocity (ft/sec) 

N = pump RPM 

oe oe flow rate per pump (fee) 

Number of pumps = 6 
WF = actual fuel flow rate per engine (lbm/hr)X100.0 


Number of engines = 6 
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AI = actual flushed inlet throat area per inlet 
(£t*)X50.0 


Number of inlets = 6 
TRST = total thrust of the six (6) waterjet system (l1bf) 
TE = engine torque (ft-lbf) 


TP 


pump torque (ft-lbf) 


Notice that the gas turbine is dually protected from 
surge, overtemperature and transient mechanical overstresses 
by 1) the fuel demand rate limiter, and 2) by the PI con- 
troller gains scheduler of Fig. 39. Thus for large errors 
the fuel demand signal received by the engine will at most 
be a six second ramp signal. This is verified by the WF 
plots in Fig. 40. 

The response in pump RPM shows that the controller 
setup of Fig. 39 is fast and stable. In addition, it has 
no overshoot. There is also no secondary oscillation due 
to the rate and hard limits of the Weg scheduler during 
eect. as this was taken care of by the negative sloping 
Piict Om Or Kon with respect to the absolute value of error. 
Thus as the error decreases the fuel pump still receives an 
increasing Weq signal. 

The flushed inlet throat area (A_) increased in the 


al 


initial transient to accommodate more flow at low S.E.S. 
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Se-ea corgenerace the acceleration force. It finally 
decreases as the S.E.S. builds up Seana eaaning War 
close to 0.8 value during transient and exactly 0.8 at steady 
state. Thus, setting a hard limit alone on the integral 
part of the PI controller for the flushed inlet is enough 
to overcome the rate and hard limit of the flushed inlet 
actuator. 

The pump torque ( T,) and engine torque (@) are 
plotted to check mechanical stress overloads during 


transients. 


(b) Deceleration from 635 to 150 RPM Throttle Setting 


This simulation was also conducted by putting a step 
change in RPM setting. In order to protect the gas turbine 
gas generator turbine inlet from sudden change in tempera- 
ture the fuel demand rate limiter for decreasing fuel flow 
rate was set arbitrarily to six second samp. The PI con- 
troller gains scheduler of Fig. 39 likewise contributes to 
this protection as an added precaution. The results are 
plotted in Fig. 41. The labels and scales are the same as 
those of the acceleration test. 

From the plots it can be seen that the response of 
pump RPM was slower than in acceleration and has slight 


overshoot. This was expected since the damping term in 





103 


GGc 


MNILLAS STLLOVHL Wdy of! OL SEI 


OST 


(SONG) 3S) 


OOt 





JSWIs 


wods 


No wYIwWIAA 


GS 


ey 


iy AyOdI4 


GOT 


ie: 


Gee 


OCYV 


Ges 


O69 


OOZ 


(Wiles iiNet 7s 


(aso 7S) ), 0 


Besea< OS) Iv 


54 





, 104 


OMILAS AlllodHlL wy ost of 569 Wodd Nollwda14a73¢ 


(SGN 0950) SAT 1 


GGe OST GOT 0S 


gly 
™“ 











dt 


LOOL 


ea!) i: i|  [ &-oe 


iy Bunoid 


Ce ee at Ce? i Sds 


Ger es) iL 


(SH/“H1*OOT) 4M 





105 


5 
oe &, 


pump RPM varies approximately to the square of pump RPM. 
Such an overshoot, however, is within acceptable limits 
and in addition it occurs at approximately 10% of the 
maximum intermittent power rating of the engine. 

Thus, it is evident from the pump RPM responses 
that the PI centroller gain of Fig. 39 works optimally 
throughout the whole operating range of the waterjet system, 
that is, 1) fast response, 2) less than 5% overshoot, 

3) stable and 4) no steady state error. 

From the plot of We it can be seen that the actual 
fuel flow is rate and hard limited even without the actuator 
itself being rate and hard limited. 

The plot of the flushed inlet throat area (A. ) showed 
an initial decrease to its minimum hard limit for several 
seconds in an attempt to bring V/V, equal to 0.8. This did 
not present a problem on the integral part of the PI con- 
troller for the flushed inlet since the integrator was also 
set to have a minimum hard limit equal to the minimum hard 
limit of the flushed inlet actuator. As the water flow rate 
Q and S.E.S. velocity ve came to a steady state, the inlet 
throat area (A. ) increased to its steady state value so that 
Wa Aae.. is exactly equal to 0.8. 

The pump torque (T.) and engine torque Ce.) were 


plotted to check back flow of power during deceleration. 
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CONCLUSIONS 


From the result of the open loop simulation to a 
step fuel flow rate for the waterjet propulsion for S.E.S. 
several conclusions can be made. The first conclusion is 
that despite the pump being designed not to cavitate at 
hump speed with a thrust to resistance ratio of 1.2, cavi- 
tation can still occur during transients since designed 
pump head and designed pump flow rates for the hump speed 
can be reached way earlier before the S.E.S. goes to hump 
speed during acceleration. However, such cavitation only 
occurs for a few seconds. 

The second conclusion is that for a fixed designed 
area nozzle the pump is operating very DOS to the designed 
Ce) and WY at steady state for all power levels; likewise, 
the gas turbine is operating very close to the point of 
maximum efficiency for all power levels at steady state. 

The third conclusion is that in the criteria for 
designing the speed at which the inlet area can change, there 
are two variables involved, namely the flow rate Q and the 
Beles. velocity Vo: During transients it was noted that Q 


enanges faster than V,. Thus, if local cavitation Ibs Tele 
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be avoided in the inlet during transient, the rate of 
change of the term WAV, should be considered. 

As already discussed in the MBODE analysis of poles 
and zeroes locations of the system the prime mover and the 
pump dynamics can be decoupled from the S.E.S. motion 
dynamics so that the design of a stable control for pump 
RPM can be simplified. A similar finding was found to be 
true also between the inlet area dynamics and the S.E.S. 
motion SoBe - 

In the fuel demand scheduled rate limiter ieee area 
in the RPM control of pump, cavitation can be avoided during 
transient by setting the rate limiter to some lower value 
for the normal operation of the S.E.S. For the hard accel- 
eration maneuvers the fuel demand scheduled rate limiter 
an be set according to what the gas turbine characteristic 
may allow for surge and overtemperature. 

The control algorithm of Fig. 39 works for a rate 
limited and hard limited actuator and is recommendable 
Since analog or digital implementation is relatively simple. 
It also works throughout the whole nonlinear operating range 


of the waterjet system. 
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RECOMMENDATIONS 


The paper presented as it is, is like a textbook 
approach to the basic principles involved in the system. 
Several engineering approximations were made for a ine 
plified digital simulation and control of the pump RPM and 
flushed inlet area. The following recommendations are 
then made for future work in this area of study. 

1. In the control of the flushed inlet, ont the 
throat area was considered in this work. This was in- 
tended to serve as a reference in the actual reshaping of 
the whole flushed inlet geometry consistent with maintaining 
maximum possible pressure recovery, avoidance of cavitation 
and minimum drag. The mechanical implementation of this 
scheme deserves some further study. 

2. The fuel demand scheduled rate limiter still 
needs some refinement to avoid cavitation during acceleration 
for S.E.S. normal operation. Such a limiter should be 
consistent with the operating height of the pump above the 
waterline, and S.E.S. motion, prime mover and pump dynamics. 

3. In order for the S.E.S. to maintain a straight 


line motion there is a need for another set of controls to 





LOg 


regulate the thrusts of all the waterjet systems. Similarly, 
turning maneuvers with and without the use of thrust de- 
flectors deserve some further study. 

4. Backing maneuvers with thrust reversers were 
not included in the study due to the absence of data on 
flushed inlet performance in the backward motion. This 
problem should also be investigated. 

5. It will be recommendable to implement a hybrid 
Simulation of the system so that the complete gas turbine 


model as presented in Chapter II can also be included. 
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APPENDIX I 


The following are the particulars of the S.E.S. 


under study: 


Length to beam ratio = 2.37 

Cushion pressure to length ratio = 2.5 
Cushion pressure = 412 ete 
Displacement = 2109 tons 


Number of lift engines = 2 


Number of propulsion engines = 6 

Designed cruise speed = 80 knots 

Maximum intermittent engine rating = 22,500 HP 
Maximum continuous engine rating = 20,500 HP 


Range = 2,000 miles 
Hump speed = 30 knots 


Bare hull resistance (smooth sea, cushion borne): 


Velocity Resistance Velocity Resistance 
(Knots) (lbf) (Knots) (lbf) 
28.45 422918 45 29257 
30 429438 50 282030 
35 367391 55 271360 
40 326606 60 265899 





Li 


PeLrocity Resistance Velocity Resistance 
(Knots) (lbf£) (Knots) (lbf) 
65 261984 35) 269772 
70 258695 90 278874 
75s) 259040 |S 290098 


80 263038 100 303252 
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APPENDIX II 


The steady operating point values of S.E.S. velocity 


(V__), pump RPM (N ), water flow rate per pump (Q06) 
oO 


OO 


Ehroat area per inlet (A; ) and redwétion gear efficiency 
fo) 


Sa ) for different values of fuel flow rate per engine 
O 


(We) are: 

ae Voo 
(lbm/hr) (ft/sec) 

800 = hee: 
2,000 PAO TA iL 
3,000 28.96 
4,000 34.37 
5,000 43.22 
6,000 46.62 
7,000 49 .67 
BOO LOS 0 
B00 124.0 
9,000 124.0 
9,300 BB oe 
10,000 142 


(RPM) 
ae 

241.8 
346.8 
410.3 
465 .6 
514 

Doom 
Shehe) eee, 
605.8 
S20-9 
5625.0 


647 .9 


2 


(£t>/sec) 


47,38 


LOdea6 


146.4 


S27 


Ot. 


2 USE 0 


23506 


2 See 


Z2OGi2 


280 .4 


29220 


298.6 


Aio 
(£7) 


62507 
6.318 
eyo ul bi) 
Geely 
Seal] 
Gmol7 
6.317 
Gsowy 
JeZO2 
Zee 2 
Zee OL 


2.628 


—_ 


0.858 


O22 

0.944 
Om a2 
OF2S7 
0.9161 
0.962 
OF IG2 
0.964 
02905 
0.966 


0.966 
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From the above operating point values, power turbine 


torque (Ge ), pump flow coefficient (D.) and head coef- 
O 


ficient (Ww ) can be determined. Thus the matrix A of the 
O 


linearized waterjet system equation (3-8) can be evaluated 


so that: 
f W = S00 
or fd 8 
-0 .0264 O20 
0.0 -0 .0293 
OES L72 1225 
> 
0.0 -4 .459 
OO O20 
0220 020 
for Weg = 5,000 
-0 .1088 0.0 
0.0 -0 .0106 
2.119 Ss ee 
A = 
0.0 -20.179 
OO 0.0 
0.20 0.0 


0.00514 


-0.00757 


-3.024 


OC 2s 9 


-0 .04479 


=12.537 


0.0 


0.0 


0.0023 02 


0.0 


-0.5 


0.0 


0.001449 


-0 .00055 
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for W__ = 10,000 
fd 


-0.03283 0.0 0.02556 0.0 
0.0 -0.1754 0.0019 0.0014 
VEL 7.704 -19.706 0.0 
0.0 -32.506 0.0 -0.5 
0.0 0.0 0.0 0.0 

O20 0.0 0.0 0.0 





The above constants were used in designing controller 


gains for small disturbance. 
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APPENDIX III 


Subroutine EQSIM is required by DYSYS for the solution 
of the dynamic equation of the system. The DYSYS program is 
available at M.I.T. mechanical and civil engineering computer 
facilities. 

Subroutine LM250 represents the LM2500 performance 
map of Fig. 8. The inputs to the subroutines are power 
turbine RPM and fuel flow rate and the output is power 


turbine torque. 
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